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*An  open-tube,  sliding-boat,  liquid-phase  system  has  been  used 
to  grow  ZnxCd-^_  -SnP2  epitaxially  on  InP  substrates.  X-ray  diffrac¬ 
tion  measurements  indicate  that  this  alloy  grows  on  ( 100}  substrates 
with  the  lattice-matched  a-axis  in  the  growth  plane  and  the  c-axis 
in  the  growth  direction.  Unintentionally-doped  layers  have  electron 
concentrations  as  high  as  3  ajJ-Q^cja-^  with  mobility  values  of  about 
2,000  cro2/V  sec.  These  mobility  values  are  substantially  larger 
than  have  been  obtained  in  the  equivalent  III-V  materials  at  simila 
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concentrations . 

Attempts  to  grow  ZnGeAs-  epitaxially  on  GaAs  substrates  with 
a  Zn-Ge-AsCl3  and  a  Zn-GeCl .-Ge-As V  vapor-phase  system  have  met  wit 
little  success.  Although  tne  first  system  gives  good  control  over 
the  As  flux  and  the  second  over  Ge,  reproducible  results  have  not 
been  achieved  in  either  system.  Our  experience  with  these  methods, 
however,  indicates  that  most  of  the  reproducibility  problems  can  be 
overcome  with  a  Zn-GeCl^-AsCl,  reactor.  We  are  currently  investi¬ 
gating  this  method. 

Recent  discussions  in  the  literature  have  pointed  out  several 
problems  with  the  conventional  theory  of  heterostructures.  We 
present  a  general,  first-order,  nonequilibrium  thermodynamic  ana¬ 
lysis  which  points  out  that  one  problem  with  conventional  theory  is 
that  it  neglects  a  generalized  thermodynamic  force.  This  force  can 
be  associated  with  the  gradient  of  the  effective  density-of-states. 
In  heterostructures  such  gradients  cause  the  charge  carriers  to 
move  into  regions  with  higher  density-of-states,  thus  increasing 
the  entropy  of  the  system  in  agreement  with  the  second  law  of  ther¬ 
modynamics.  We  also  show  that  if  the  electrochemical  potential  is 
referenced  to  the  infinite  vacuum  level,  then  the  electron  affinity 
rule  must  obviously  be  valid. 

^  Calculations  are  presented  on  the  I-V  characteristics  of  a 
space-charge  limited  diode.  If  the  electrons  exhibit  inertial  motio 
(are  not  scattered) ,  then  they  will  quickly  reach  nonparabolic 
regions  of  the  conduction  band.  When  this  nonparabolicity  is  taken 
into  account,  the  current  appears  ohmic  rather  than  exhibiting  a 
V3/2  behavior,  as  has  been  previously  supposed. 
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II-IV-V2  CHALCOPYRITES  FOR  HIGH-SPEED  DEVICES 

1.  RESEARCH  OBJECTIVES 

Because  of  their  lower  effective  masses,  different 
band  structure,  and  good  lattice  match  to  their  analog 
III— V  compounds,  some  of  the  II-IV-v2  compounds  may  be 
important  for  heterostructures  and  high-speed  device  ap¬ 
plications.  The  objectives  of  the  work  reported  here  are 
to: 

(1)  develop  techniques  for  the  growth  of  high-quality 
epitaxial  Zn  Cd.  SnP~  on  InP  and  GaAs  and  ZnGeAs7 
on  GaAs  or  Ge; 

(2)  minimize  electrically-active  defects  and  residual 
impurities  in  these  materials  to  achieve  well- 
controlled  doping; 

(3)  investigate  the  electrical  transport  properties 
of  these  materials  to  determine  their  suitablilty 
for  heterostructures  and  high-speed  device  appli¬ 
cations;  and 


(4)  develop  the  technology  required  for  the  fab¬ 
rication  of  suitable  device  structures. 

With  these  results  we  expect  to  gain  substantial  insight 
into  the  potential  of  these  II-IV-V2  chalcopyrites  for 
electronic  applications. 


2.  LIQUID  PHASE  EPITAXIAL  GROWTH  OF  ZnxCd1_xSnP2 


CdSnP2  is  a  potentially  useful  compound  semiconductor 
because  it  has  about  half  (0.035m)  the  electron  effective 
mass  of  its  III-IV  analog  InP  with  a  comparable  bandgap 
(1.17  eV)  [l,2].  Previous  attempts  to  grow  CdSnP2  epi¬ 
taxially  on  InP  by  LPE  have  met  with  limited  success  due 
in  part  to  the  0.5%  a-to-a  mismatch  between  the  two 
materials  [3,4].  The  purpose  of  the  work  in  this  section 
is  to  develop  the  technology  to  produce  device  quality 
epitaxial  layers  to  see  if  the  low  electron  effective  mass 
of  this  material  can  be  utilized. 

2.1  MATERIAL  CONSIDERATIONS 

A  cubic  unit  cell  of  the  CdSnP2  chalcopyrite  structure 
is  compared  to  two  cubic  unit  cells  of  the  InP  sphalerite 
structure  in  Figure  1.  These  crystal  structures  are 
similar  with  the  column  III  element  sublattice  of  the 
sphalerite  structure  being  replaced  by  column  II  element 
and  column  IV  element  sublattices  in  the  chalcopyrite 
structure.  Because  of  the  structure  of  these  sublattices 
the  volume  of  the  chalcopyrite  cubic  unit  cell  is  approxi¬ 
mately  twice  that  of  the  sphalerite  cubic  unit  cell.  The 
chalcopyrites  exhibit  uniaxial  compression  so  that  two 
lattice  constants  are  needed  to  describe  the  lattice:  a 
and  c,  where  c/a  <_  2  for  all  chalcopyrites.  CdSnP2  is  an 
analog  of  InP  since  Cd  and  Sn  straddle  In  on  the  periodic 


T 
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table.  Some  properties  of  InP,  CdSnP2,  and  ZnSnP2  are 
presented  in  Table  1. 

The  0. 5%  a- to-a  mismatch  of  CdSnP2  on  InP  has  proven 
to  be  a  serious  impediment  to  the  growth  of  device  quality 
epitaxial  layers.  A  typical  CdSnP2  layer  grown  on  InP 
produced  in  the  early  stages  of  this  work  is  shown  in 
Figure  2.  The  surface  of  the  grown  layer  is  exceedingly 
rough  and  the  polycrystalline  nature  of  the  layer  is  re¬ 
vealed  in  the  cross-sectional  view.  The  thin  connected 
layer  at  the  interface  is  a  layer  of  regrown  InP. 

To  reduce  the  lattice  mismatch  in  this  heterostructure, 
we  decided  to  grow  the  lattice  matched  alloy  of  ZnxCd1_xSnP2 
although  other  chalcopyrite  alloys  could  be  used,  ZnSnP2 
was  chosen  because  the  lattice  matched  composition  in  this 
case  is  closer  to  CdSnP2  than  in  other  possible  alloys. 
Because  of  the  uniaxial  compression  of  the  chalcopyrites 
it  is  not  possible  to  lattice  match  both  the  a  and  c 
lattice  constants  simultaneously.  Thus,  growth  is  re¬ 
stricted  to  {100}  oriented  InP  substrates  so  that  the  a 
lattice  constants  can  be  matched  at  the  interface  with  the 
c  lattice  constant  extending  normal  to  the  growth  interface. 
Such  growth  is  expected  from  x-ray  diffraction  measurements 
of  CdSnP2  grown  on  {100}  InP.  Assuming  Vegard's  law,  the 
calculated  alloy  composition  to  lattice  match  InP  is  Zn_  . .. 

U  t  X 

Cd0. 87  SnP2  * 
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2.2  GROWTH  PROCEDURE 

The  system  shown  schematically  in  Figure  3  was  used 
for  the  epitaxial  growth.  This  system  employs  a  standard 
LPE  graphite  boat  modified  so  that  the  substrate  sits  in 
the  base  with  the  melt  wells  extending  through  the  slider. 
This  boat  was  designed  so  that  growth  experiments  could  be 
carried  out  despite  the  lack  of  phase  diagram  information 
for  this  system.  It  also  allows  us  to  initiate  growth  with 
a  saturated  solution  despite  the  loss  of  volatile  consti¬ 
tuents  from  the  melt  during  homogenization. 

During  homogenization  the  slider  is  positioned  so  that 
its  melt  wells  are  aligned  with  the  holes  in  the  top 
section.  In  this  position  the  substrate  is  covered  by  the 
slider  to  reduce  the  amount  of  volatile  constituents 
reaching  its  surface.  These  contitutents  react  with  the 
substrate  causing  severe  pitting.  Excess  source  materials, 
which  have  been  added  to  the  melt,  float  to  the  surface 
to  keep  the  solution  saturated  as  volatile  constituents 
evaporate.  Growth  is  initiated  by  separating  the  melt  con¬ 
tained  in  the  slider  from  the  rest  of  the  solution.  This 
effectively  removes  the  excess  source  material  before  the 
melt  is  positioned  over  the  substrate.  In  this  position  the 
growth  melt  is  covered  by  the  top  section  to  reduce  deple¬ 
tion  due  to  evaporation. 

The  growth  parameters  used  in  the  system  for  epitaxial 
growth  are  presented  in  Table  2.  The  source  materials  used 


Figure  3  System  used  for  LPE  grow 
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Table  2 


Typical  growth  parameters  for  the  Zn-Cd-SnP,-Sn 
system. 


SOURCE  MATERIALS:  Zn,  Cd,  SnP3,  Sn 


Melt  Composition 


CdSnP2 

19.85 

Molar 

ZnSnP2 

0.25 

Molar 

Sn 

77.50 

Molar 

P 

2.40 

Molar 

ZnSnP 

2/CdSnP2 

*  0.0125 

Substrates:  {100}  Fe  Doped  InP 


Homogenization  Temperature: 


540°C 


Cooling  Rate: 
Growth  Duration: 
in  aitu  Etch: 
Growth  Atmosphere: 
Flow  Rate : 


5° C/hour 
10  to  22  hours 
Pure  Sn 
UHP  H2 
30  ml/min 
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are  5Ns  purity  Zn,  Cd,  and  Sn  with  SnP3  used  as  a  source 
of  phosphorus.  The  SnP^  was  synthesized  in  a  sealed  tube 
system  using  5N  Sn  and  red  P^.  The  melts  had  an  initial 
equivalent  composition  of  19.85%  CdSnP2/  0.25%  ZnSnP2, 

2.40%  excess  P,  and  77.50%  Sn.  This  gives  a  ZnSnP2  to 
CdSnP2  ratio  in  the  melt  of  0.0125  which  is  about  a  tenth 
of  the  ratio  required  in  the  solid  to  lattice  match  InP. 

A  homogenization  and  initial  growth  temperature  of 
540°C  was  used.  For  an  initial  growth  temperature  of  530°C 
the  growth  was  often  polycrystalline.  It  should  be  noted 
that  these  temperatures  are  quite  close  to  the  570°C  maxi¬ 
mum  melting  point  of  CdSnP2-  After  the  growth  melt  is  posi¬ 
tioned  over  the  substrate,  the  furnace  is  cooled  at  a  linear 
rate  of  5°C/hour.  cooling  rates  of  10°C/hour  produced  poly¬ 
crystalline  layers  and  rates  of  4°C/hour  produced  no  growth 
due  to  melt  depletion  by  evaporation.  Cooling  times  be¬ 
tween  10  and  22  hours  were  used.  A  10  hour  cooling  time 
'  gave  lym  to  2ym  of  growth,  while  a  22  hour  time  gave  around 
15ym  of  growth. 

In  most  instances  an  in  situ  etch  of  the  substrate  in 
pure  Sn  was  made  prior  to  growth.  A  typical  epitaxial 
layer  grown  using  these  parameters  is  shown  in  Figure  4. 

As  can  be  seen  from  this  figure  Znxc<^i_xSnP2  ePitaxial 
layers  are  much  smoother  than  the  CdSnP2  layers  grown  on 
InP,  although  they  are  still  not  of  device  quality.  The 
black  regions  in  the  view  of  the  upper  surface  are  etch 
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pits  which  show  up  as  notches  in  the  cross-sectional  view. 
These  pits  extend  to  various  depths  into  the  epitaxial 
layers  and  often  extend  to  the  interface. 

The  sharp  dependence  of  the  layer  morphology  on  zinc 
concentration  in  the  melt  is  shown  in  Figure  5.  This  figure 
shows  two  epitaxial  layers  grown  in  the  same  manner  except 
for  the  zinc  concentration  in  the  melts.  The  layer  shown 
in  Figure  5a  was  grown  from  a  standard  melt  and  that  shown 
in  Figure  5b  was  grown  from  a  melt  having  a  6%  larger  mass 
of  zinc.  This  latter  layer  is  polycrystalline. 

2.3  ELECTRICAL  MEASUREMENTS 

Van  der  Pauw  measurements  were  obtained  for  the  epi¬ 
taxial  layers  with  the  results  shown  in  Figure  6.  All 

these  layers  were  n  type.  Most  had  room  temperature  cal- 

18  —  3 

culated  carrier  concentrations  around  10  cm  and  mobil- 

2 

ities  around  100  cm  /V  sec,  similar  to  layers  of  CdSnP2 

-2 

grown  on  InP.  The  resistivities  were  around  3  x  10  ft  cm. 
Notice  there  are  no  substantial  changes  in  these  values  for 

measurements  made  at  77  K.  One  layer  had  a  mobility  of 

2  19-3 

2000  cm  /V  sec  with  a  carrier  concentration  of  3  x  10  cm 

-4 

The  resistivity  of  this  sample  was  1  x  10  fi  cm.  This 
sample  is  of  interest  because  the  mobility  is  about  twice 
that  which  can  be  obtained  for  equivalent  III-V  compounds 
at  similar  carrier  concentrations.  This  sample  indicates 
one  potential  advantage  of  the  low  effective  mass  II-IV-V2 
chalcopyrites  over  the  III-V  sphalerite  semiconductors. 


Figure  5  Epitaxial  layers  gr< 
composition,  and  (b 
6%  greater  Zn  mass. 


o  300  K 
v  77  K 
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Figure  6  Results  of  Hall  measurements  for  Zn  Cd,  SnP7  epitaxial  layers 
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That  is,  higher  mobilities  can  be  achieved  at  similar 
carrier  concentrations. 

The  solid  line  in  Figure  6  shows  the  theoretical  de¬ 
pendence  of  the  mobility  on  the  carrier  concentration  for 
scattering  from  ionized  impurities  with  screening.  This 
curve  was  calculated  from  the  Brooks-Herring  model  [5] 
using  degenerate  statistics  to  allow  for  the  high  carrier 
concentrations.  The  concentration  of  ionized  impurities, 

Nj,  in  this  model  was  used  as  the  only  adjustable  parameter. 

20  -3 

To  fit  the  high  mobility  data,  a  value  of  Nj  =  1  x  10  cm 
was  used.  As  can  be  seen,  there  is  reasonable  agreement 
between  the  experimental  data  and  the  theoretical  curve. 

This  indicates  that  the  decrease  in  mobility  with  decreas¬ 
ing  carrier  concentration  is  mostly  due  to  heavy  compensation . 

Because  high  conductivity  in  homogenieties  in  an  epi¬ 
taxial  layer  can  give  anomalously  high  calculated  mobilities 
[6] ,  further  experiments  were  performed  on  the  high  mobility 
sample  to  determine  whether  the  mobility  value  was  real. 
First  this  sample  was  cleaved  into  four  pieces  and  the 
pieces  remeasured.  The  effects  of  a  conducting  inhomo- 
geniety  on  resistivity  and  Hall  measurements  can  usually  be 
eliminated  in  this  manner.  The  mobilities  and  carrier  con¬ 
centrations  of  each  piece,  however,  were  very  close  to  those 
of  the  original  sample.  This  can  be  seen  from  the  tight 
grouping  of  the  data  in  the  high  concentration  region  of 
Figure  6.  To  further  substantiate  these  data,  measurements 


of  the  Hall  constant  as  a  function  of  magnetic  field  were 
performed-  These  results  are  plotted  in  Figure  7,  where 
a  decrease  in  the  Hall  constant  with  increasing  magnetic 
field  is  shown.  This  behavior  is  typical  of  homogeneous 
n-type  material.  If  there  had  been  a  conducting  inhomo- 
geniety,  such  as  those  which  give  anomalously  high  calcu¬ 
lated  mobilities,  the  Hall  constant  would  have  increased 
with  increasing  magnetic  field.  Thus,  the  mobility  and 
carrier  concentration  of  this  epitaxial  layer  are  apparently 
not  due  to  inhomogeniety  effects  and  appear  to  be  an  ac¬ 
curate  indication  of  the  high  mobility  which  can  be  attained 
in  this  material. 

The  high  ionized  impurity  concentration  in  these 
samples  may  be  due  to  loss  of  volatile  constituents  from 
the  epitaxial  layers  during  growth.  This  is  supported  by 
the  results  of  high  temperature  Hall  measurements  on  SiC>2 
encapsulated  platelets  of  CdSnP2.  These  measurements  showed 
a  decreasing  mobility  and  increasing  carrier  concentration 
with  annealing  time.  Annealing  of  the  layers  in  a  phos¬ 
phorus  or  phosphorus  and  cadmium  overpressure  should  then 
give  increasing  mobilities  and  decreasing  carrier  concen¬ 
trations  as  Nj  is  reduced.  Such  experiments  are  the  subject 
of  a  future  investigation. 

2.4  CONCLUSIONS 

The  method  presented  for  the  LPE  growth  of  the  lattice 
matched  alloy  Zn  Cd.  SnP-  on  InP  has  yielded  epitaxial 

X  X  ^ 
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layers  substantially  better  than  layers  of  CdSnP2  on  InP. 
The  morphology  of  these  layers.  However,  is  not  of  device 
quality.  The  reasons  for  this  poor  morphology  may  be: 

(1)  zinc  depletion  from  the  melt  due  to  its  high  segre¬ 
gation  coefficient;  (2)  constitutional  supercooling;  and/or 
(3)  antiphase  boundaries  due  to  the  additional  sublattices 
of  the  layer  compared  to  the  substrate.  We  are  currently 
investigating  these  possibilities. 

Van  der  Pauw  measurements  on  the  epitaxial  layers  have 
shown  that  the  mobility  of  this  material  appears  to  be  sub¬ 
stantially  higher  than  can  be  obtained  in  equivalent  III-V 
compounds  at  similar  carrier  concentrations.  Further  work 
needs  to  be  done  to  substantiate  this  finding.  We  plan  to 
perform  additional  growth  experiments  as  well  as  annealing 
experiments  to  determine  if  these  higher  mobility  values 
can  be  obtained  at  lower  carrier  and  ionized  impurity 
concentrations . 
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3.  VAPOR  PHASE  EPITAXIAL  GROWTH  OF  ZnGeAS2 

The  tetragonal  chalcopyrite  compound  ZnGeAs2  with  a 
bandgap  of  1.15eV  is  a  structural  and  electronic  analog 
of  GaAs.  If  it  could  be  reproducibly  prepared  in  single 
crystalline  form  and  suitably  doped,  ZnGeAs.,  could  be  an 
interesting  material  for  microwave  devices .  In  this 
section  we  discuss  vapor  phase  techniques  which  have  been 
used  in  an  attempt  to  achieve  this  purpose. 

3 . 1  GROWTH  PROCEDURE 

Attempts  to  synthesize  ZnGeAs2  were  carried  out  by 
employing  high  purity  materials  Zn,  Ge,  As,  and  GeCl4< 

The  substrates  were  semi-insulating  Cr-doped  GaAs. 

Figure  8  shows  the  reactor  which  was  used  for  this  purpose. 
It  consists  of  a  silica  tube  in  a  furnace  with  three  in¬ 
dependently  controllable  zones.  The  left  hand  zone  con¬ 
trols  the  As  vapor  pressure,  the  center  zone  controls  the 
Zn  vapor  pressure,  and  the  right  hand  zone  is  the  reaction 
zone.  The  GeCl^  is  contained  in  a  bubbler  held  at  0  to 
10°  C  and  H2  is  used  as  the  carrier  gas.  An  additional 
H2  flow  is  also  provided  to  flush  the  reaction  tube  and  to 
sweep  the  reactants  into  the  reaction  zone . 

The  GaAs  wafers  used  in  the  experiments  were  lapped 
smooth  and  etched  vigorously  in  an  etch  consisting  of 
5H2S04 : 1H202 :1H20  to  produce  a  bright,  mirror  finish. 
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The  H2  flow  rate  for  GeCl^  was  100  ml/min.  The  additional 
H2  flow  rate  was  100  ml/min. 

To  try  to  find  the  experimental  conditions  under 
which  ZnGeAs2  could  be  grown,  the  Zn  and  As  temperatures 
were  set  at  600  and  580°  C,  respectively.  These  tempera¬ 
tures  provided  a  large  enough  flux  of  these  species  to 
essentially  flood  the  reactor  zone.  With  the  Ge  tempera¬ 
ture  set  at  560°  C  to  produce  the  disproportionation 
reaction  with  GeCl4,  the  temperatures  of  the  GeCl4  and  the 
reaction  zone  were  systematically  varied.  Over  most  of  the 
range  of  these  experiments,  little  or  no  growth  occurred 
on  the  {100}  GaAs  substrates.  With  substrate  temperatures 
around  700°  C  and  GeCl4  temperatures  in  the  range  from  4  to 
10°  C,  a  number  of  good  quality  epitaxial  layers  were 
obtained.  Subsequent  analysis,  however,  indicated  that 
these  layers  were  mostly  Ge. 

3.2  EXPERIMENTAL  RESULTS 

Figure  9  shows  the  upper  surface  and  interface  of 
one  of  these  Ge-GaAs  heterostructures.  Resistivity  and 
Hall  measurements  on  these  Ge  epitaxial  layers  on  semi- 
insulating  GaAs  substrates  indicated  very  high  n-type 

conductivity.  Typical  values  of  electron  concentration 

19-3  2 

and  mobility  were  5  x  10  cm  and  100  cm  /V  sec,  re¬ 
spectively.  Some  of  these  n+  Ge  layers  were  also  grown 
15-3 

on  10  cm  n-type  GaAs.  The  specific  contact  resistance 
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(b)  Interface 


Figure  9  Ge-GaAs  heterostructure  grown  in  the  GoCI,  system 
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2 

of  these  n-n+  heterostructures  was  about  0 . lftcm  .  Although 
these  Ge  layers  are  quite  good,  we  were  not  able  to  grow 
the  compound  ZnGeAs2  under  these  conditions. 

Since  the  preparation  of  the  compound  appears  to  be 
so  critically  dependent  upon  the  reaction  conditions, 
x-ray  diffraction  measurements  were  made  to  determine  the 
composition  of  the  material  deposited  at  various  places  in 
the  reactor.  Several  growth  runs  were  made  with  different 
growth  parameters.  The  materials  deposited  out  at  the 
different  places  of  the  furnace  tube  were  prepared  for 
diffraction  analysis  by  grinding  them  into  a  fine  powder 
and  placing  them  on  a  flat  plate.  This  plate  was  then 
inserted  into  the  specimen  holder  of  the  diffractometer. 
Several  diffraction  patterns  were  obtained.  The  materials 
were  then  identified  from  their  diffraction  data  with  the 
aid  of  a  Hanawalt  index  and  file.  Some  of  these  results 
are  shown  in  Table  3. 

Typically,  the  factors  which  affected  the  nature  of 
the  deposits  most  were  substrate  temperature  and  reactant 
ratios.  Polycrystalline  As  was  formed  downstream  at  lower 
temperatures  while  Ge  and  Zn^AS2  were  formed  at  higher 
temperatures.  The  Ge/I^  ratio  or  Ge  input  flux  variations 
were  achieved  by  changing  the  vapor  pressure  at  the  GeCl4 
source  while  keeping  the  total  flow  constant.  The  data 
show  a  strong  Ge  growth  rate  dependence  with  Ge/I^  ratio. 
The  choice  of  reactant  vapor  pressure  ratio,  Zn :GeCl4 : As4 , 
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appears  to  be  important.  When  the  ratio  is  1:1:2,  large 
quantities  of  Zn3As2  platelets  form  together  with  Ge 
platelets . 

3.3  CONCLUSIONS 

The  preparation  of  ZnGeAs2  with  the  Zn-GeCl4-Ge-As4 
system  was  found  to  be  much  too  strongly  dependent  upon 
the  various  growth  parameters.  When  As4  was  transported 
over  the  Zn  boat,  Zn3As2  platelets  were  formed.  Apparent¬ 
ly,  this  greatly  reduces  further  transport  of  Zn  or  As 
vapor  and  results  in  the  depletion  of  Zn  or  As  in  the  reac¬ 
tion  zone.  This  appears  to  be  the  main  problem  with  this 
Zn-GeCl4~Ge-As4  system:  it  gives  excellent  control  of 
the  Ge  flux,  but  very  little  control  for  Zn  and  As. 

Previous  to  this  study,  we  investigated  a  Zn-Ge-AsCl3 
reactor  system  for  the  synthesis  of  ZnGeAs2  and  found 
that  it  produced  good  As  control.  For  these  reasons  we 
are  currently  attempting  to  synthesize  the  compound  in 
a  Zn-GeCl4-AsCl3  reactor,  which  incorporates  the  best 
features  of  the  two  previous  systems  without  their 
disadvantages . 

This  new  system  is  shown  in  Figure  10.  The  reactor 
consists  of  a  44  mm  (OD)  quartz  tube  having  separate 
reactant  inlets  and  a  furnace  with  three  independently 
controllable  zones.  The  flow  rates  of  Zn,  GeCl4 ,  and 
AsCl3  are  adjusted  by  independent  carrier  gas  flowmeters. 


The  left  hand  zone  controls  the  Zn  vapor  pressure,  while 
the  right  hand  zone  is  the  reaction  zone.  A  baffle  in 
the  Zn  chamber  is  used  to  isolate  the  Zn  from  the  As. 

The  AsCl3  and  GeCl4  are  kept  in  external  bubblers.  The 
primary  advantages  of  this  system  are:  (1)  it  ensures 
that  the  constituents  do  not  react  until  they  reach  the 
reaction  zone;  (2)  it  provides  better  over  all  control 
of  the  three  constituents;  and  (3)  it  utilizes  only  one 
solid  source  which  makes  it  much  simpler  to  operate  than 
the  previous  reactor  systems. 


4.  NONEQUILIBRIUM  THERMODYNAMICS  OF  HETEROJUNCTION 


The  recent  interest  in  heterostructures  for  lasers, 
transistor  emitters,  real-space  transfer  devices,  quantum 
well  lasers,  modulation-doped  superlattices,  and  other  de¬ 
vices,  has  prompted  a  re-examination  of  basic  heterojunction 
theory.  To  understand  quantitatively  the  operation  of 
these  devices  it  is  necessary  to  determine  how  the  discon¬ 
tinuity  in  the  energy  gaps  of  a  heterostructure  is  divided 
between  the  conduction  and  valence  bands : 

Aeg  =  Aec  +  Aev  (1) 

The  basic  assumption  of  heterostructure  theory  is  the  "elec¬ 
tron  affinity  rule":  that  is,  the  conduction  band  discon¬ 
tinuity  is  simply  the  difference  between  the  electron 
affinity  of  the  two  materials, 

Aec  =  *i  "  X2  /  (2) 

presumably  referenced  to  the  infinite  vacuum  level. 

4 . 1  CURRENT  PROBLEMS 

Kroemer  [7]  has  pointed  out  a  number  of  problems  with 
this  concept,  the  most  important  of  which  are:  (1)  assum¬ 
ing  the  rule  is  valid,  the  band  discontinuities  need  to  be 
known  to  an  accuracy  of  about  a  kT,  while  the  affinities, 
usually  referenced  to  a  near  vacuum  level,  cannot  be  meas¬ 
ured  that  accurately;  (2)  the  basic  differences  between 
free  surfaces  and  solid  interfaces  are  so  great  that  there 
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is  no  fundamental  reason  for  expecting  any  relationship 
between  electron  affinities  and  band  discontinuities.  Thus, 
to  resolve  this  problem,  it  is  necessary  to  forget  about 
phenomenological  surface  properties  and  look  at  the  inter¬ 
face  problem  from  a  microscopic  point-of-view. 

In  spite  of  these  arguments,  Shay  et  al  [8]  compared 
measured  band  discontinuities  for  several  hetero junctions 
with  measured  electron  affinities  from  photoemission  data 
and  calculated  affinities  from  the  dielectric  two-band  model. 
They  found  that  the  agreement  between  the  measured  discon¬ 
tinuities  and  those  calculated  from  the  affinity  rule  was 
very  good.  In  addition,  Frensley  and  Kroemer  [9]  found  from 
a  pseudopotential  calculation  of  a  hetero junction  interface 
that  the  affinity  rule  is  approximately  correct.  Although 
this  is  a  somewhat  surprising  result,  Kroemer' s  original 
objections  are,  of  course,  still  valid. 

More  recently,  Adams  and  Nussbaum  [10]  have  examined 
the  affinity  rule  from  an  electrostatic  point-of-view.  They 
suggest  that  instead  of  keeping  the  vacuum  level  continuous 
(the  electron  affinity  rule)  it  may  be  more  appropriate  to 
assume  continuity  of  the  intrinsic  level.  This  paper  is 
interesting  because  it  also,  although  unintentionally,  dem¬ 
onstrates  another  problem  with  hetero junction  theory:  the 
usual  equilibrium  and  non-equilibrium  formulations  in  the 
depletion  approximation  are  not  self-consistent.  This 
inconsistency  occurs  regardless  of  which  basic  assumption 


is  used. 


The  problem  arises  in  the  following  manner:  For  sim¬ 
plicity  assume  an  abrupt  p-n  hetero junction.  In  equilibrium, 
the  built-in  potential  of  the  space-charge  region  is  given 
by, 


where  and  Na  are  the  donor  and  acceptor  concentrations  on 

the  n  and  p  side  of  the  hetero junction,  respectively,  and 

n.  and  n.  are  the  intrinsic  concentrations.  When  the 
in  ip 

hetero junction  is  forward-biased  with  a  voltage  V,  in  the 
depletion  approximation  the  holes  injected  into  the  n-side 
at  the  edge  of  the  depletion  region  is  usually  found  to  be, 

r  g^-vi 

Pn  =  PpexP[ - i?F-J  ‘  (4) 

If  V  is  turned  off,  the  excess  holes  disappear  and 


2 
n . 


Comparing  Eqs .  (3)  and  (5),  we  see  that  this  analysis  is 
only  valid  for  a  p-n  homo junction. 

The  problem  we  wish  to  consider  is  the  latter  one, 
although  the  electron  affinity  rule  is  a  natural  consequence 
of  our  analysis.  We  will  show  that  the  inconsistency  arises 
from  the  neglect  of  a  thermodynamic  force  term,  which  is  well 
understood  but,  apparently,  not  widely  known. 
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4.2  THERMODYNAMIC  ANALYSIS 

For  any  solid  it  can  be  shown  [11]  that  the  entropy 
source  strength  (entropy  production  rate  per  unit  volume)  is, 


a  =  J  •  V 
s  Q 


In  this  equation  Jq  is  the  heat  flux  (flux  of  internal  ener- 
qv) ,  j ^  are  the  particle  fluxes  for  each  i  component,  are 
the  forces,  ik  are  the  chemical  potentials,  A^  are  the  chemi¬ 
cal  affinities  of  the  j  quasi-chemical  reactions,  and  are 
the  reaction  rates.  From  the  second  law  of  thermodynamics, 
the  entropy  source  strength  is  positive  definite  and  equal 
to  zero  in  equilibrium. 

The  basic  assumption  underlying  this  model  is  the  con¬ 
cept  of  "local  equilibrium".  That  is,  although  the  whole 
system  may  be  far  from  equilibrium,  the  particles  in  each 
differential  volume  are  in  equilibrium  with  each  other.  If 
one  considers  the  small  distances  and  interaction  times  among 
particles  within  a  differential  volume  and  the  large  dis¬ 
tances  and  interaction  times  among  particles  from  different 
differential  volumes  of  a  system,  the  concept  of  local  equi¬ 
librium  is  a  reasonable  one. 

For  our  purposes  here  we  neglect  any  temperature  gradi¬ 
ents  or  quasi-chemical  reactions  and  assume  that  the  only 
force  on  the  particles  is  an  electrostatic  potential  gradi¬ 
ent.  Equation  (6)  then  reduces  to. 
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°s  =  '  o,  <?> 

i 

where  is  the  electrochemical  potential.  Consider  the 
form  of  as  near  equilibrium  by  expanding  it  in  a 
series  of  the  generalized  forces  around  V?i  =  0.  The  first 
term  is  zero  because  og  =  0  in  equilibrium,  and  the  second 
term,  linear  in  V£^,  is  zero  because  og  is  positive  definite. 
The  third  term,  quadratic  in  V^,  is  the  first  non-zero  term. 
Retaining  only  this  term  we  can  see  that,  to  first  order, 
near  equilibrium. 


a  - 


(8) 


Going  from  particle  current  to  electric  current,  Eq. 
(8)  tells  us  that  the  correct  expression  for  the 
ecruation-of-state  of  this  isothermal  system  is. 


J.  = 

i 


n.  u • . 

ii  i 


(9) 


where  n^  and  vk  are  the  concentration  and  mobility  of  the  i 
particles.  Equation  (9),  of  course,  is  a  well-known  result. 
We  went  through  this  analysis  to  show  that  it  is  correct  to 
first  order,  under  the  assumptions  discussed,  and  that  it 
can  be  derived  from  the  second  law  of  thermodynamics.  In 
general,  all  of  the  quantities  in  Eq.  (9)  can  vary  with 
position  and  time. 


I 


We  now  need  to  obtain  an  expression  for  the  electro¬ 
chemical  potential  for  indistinguishable  particles.  For  a 


-34- 


grand  canonical  ensemble  of  particles  near  equilibrium,  it 
can  be  shown  [12]  that  the  distribution  function  is 


approximately , 


f  =  1  +  exp 


(It)] 


That  is,  it  has  the  same  form  as  the  Fermi -Dirac  distribution 
except  that,  in  general,  f,  e,  t. ,  and  T  can  vary  in  time  and 
position.  Equation  (10)  is  obtained  under  the  local  equili¬ 
brium  assumption  and  is  valid  only  to  first  order.  Thus, 

Eqs.  (9)  and  (10)  are  derived  with  the  same  assumptions  and 
approximations,  and  we  can  use  them  to  analyze  heterostruc¬ 
tures  keeping  this  in  mind. 

4.3  HETOROJUNCTION  MODEL 

Let  us  apply  Eq.  (10)  to  obtain  an  expression  between 
the  concentration  of  electrons  in  the  conduction  band  and 
their  electrochemical  potential.  In  the  usual  way,  we 
multiply  the  distribution  function  times  the  density-of-states 
and  integrate  over  the  band  to  obtain. 


g(e)f(e)de  , 


where  e  is  the  bottom  of  the  conduction  band.  For  a  dilute 
c 

distribution  of  electrons,  Eq.  (11)  gives, 


n  =  Nc  exp  (£  -  ec) , 


where  N  is  the  effective  density-of-states  of  the  conduc- 
c 

tion  band.  In  general,  all  four  quantities  in  Eq.  (12)  can 


vary  with  time  and  position. 


-  ** 
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Using  Eqs .  (9)  and  (12),  we  have  the  general 
near-equilibrium  equation-of-state  for  electron  flow  in  a 
hetero junction , 


n 


=  ny 


n’(Ec  ♦  XT  In 


(13) 


Notice,  that  all  of  the  quantities  in  this  equation  can,  in 

general,  vary  with  position  and  time  including  N  .  It  is 

*► 

the  neglect  of  the  term  -V ( lnNc)  which  leads  to  the  incon¬ 
sistency  mentioned  previously.  A  similar  expression  can  be 
obtained  for  the  hole  current: 


3P  ■  %'(£v  - kTln^)  •  (14) 

Let's  examine  each  of  the  thermodynamic  forces  in  Eqs. 

-*■  -*■ 

(13)  and  (14) ,  The  forces  7e  and  Ve  represent  variation 

v»  V 

in  composition  or  bandgap  as  well  as  applied  and  built-in 

electric  fields  due  to  space  charge.  The  component  of  these 

forces  which  has  equal  and  opposite  effect  on  electrons  and 

holes  can  be  associated  with  electric  field,  with  the  rest 

■+ 

associated  with  band-gap  variations.  The  forces  V  (Inn)  and 
-► 

-V  Xlnp)  are  the  usual  diffusion  terms.  The  thermodynamic 
forces  -V(lnNc)  and  V  { lnN^.)  can  be  regarded  as  entropy 
forces:  that  is,  they  act  to  drive  the  electrons  and  holes 
to  regions  of  higher  densi ty-of-state  or  higher  entropy. 
These  forces  are  required, therefore ,  by  the  second  law  of 
thermodynamics,  and  must  be  included  in  a  heterostructure 
analysis . 
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We  will  now  show  that  including  these  entropy  forces 
eliminates  the  inconsistancy  in  the  analysis  of  an  abrupt 
p-n  heterojunction.  In  equilibrium,  from  Eqs.  (13)  and  (14) 
the  built  in  potential  of  the  space  charge  region  is. 


(15) 


(16) 


Using  Eq.  (1)  these  two  expressions  can  be  shown  to  be  equal. 
When  the  junction  is  forward-biased,  the  holes  injected  into 
the  n-side  in  the  depletion  approximation  are. 


=  P 


N_ 
P  N 


vn 


exp 


vp 


q(V-VQ)-Ae^ 

3cT 


(17) 


Turning  the  applied  voltage  off  in  Eq.  (17) ,  we  obtain  Eq. 

(16)  and  the  inconsistancy  is  removed. 

Notice  that  in  this  analysis  we  have  made  no  assumptions 
about  the  continuity  of  the  infinite  vacuum  level  (electron 
affinity  rule)  or  any  other  level.  However,  if  we  refer  to 
Eqs.  (7)  and  (8),  in  equilibrium  there  is  no  entropy  produc¬ 
tion  and 


E(vq)2  =  0 


(18) 


or 


0 


(19) 
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If  the  electrochemical  potentials  are  referenced  to  the 
infinite  vacuum  level/  then  it  must  be  continuous. 
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5.  INERTIAL  TRANSPORT  IN  NONPARABOLIC  BANDS 

Since  the  primitive  unit  cell  of  the  II-IV-V2  chalcopy- 
rites  is  four  times  as  large  as  that  of  their  III-V  sphale¬ 
rite  analogs,  the  first  Brillouin  zone  is  one-fourth  as 
large.  For  this  reason  and  because  of  their  small  effective 
masses,  a  potential  application  for  these  chalcopyrites  may 
be  in  inertial  transport  devices  (if  such  devices  can  be 
realized) . 

The  reprint  in  this  section  shows  what  current-voltage 
characteristics  a  space-charge  limited  diode  would  have  if 
the  charge  carriers  are  not  scattered  in  a  nonparabolic  band. 
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The  expected  current  is  reduced  significantly  below  the  constant  ef¬ 
fective  mass  value  for  voltages  in  excess  of  the  bandgap,  and  the  high- 
voltage  asymptotic  relation  is  I'  rather  than 


In  devices  that  use  high-mobility,  low  effective  mass  semi¬ 
conductors,  the  mean  free  time  between  collisions  of  a  carrier 
can  become  equal  to,  or  longer  than  an  oscillation  period,  or 
the  mean  free  path  between  collisions  equal  to  or  longer  than 
a  typical  device  dimension. 

Under  such  conditions  the  carrier  transport  in  semicon¬ 
ductors  is  not  governed  by  the  traditional  friction-dominated 
Ohm’s  law  behavior.  The  motion  of  the  carriers  for  the  short 
times  or  distances  of  primary  interest  for  device  operation  may 
be  more  accurately  described  by  the  inertial  or  ballistic  equa¬ 
tion  of  motion. 

The  simplest  device  for  studying  inertial  transport  is  the 
short  space-charge  limited  semiconductor  diode.  Calculations 
of  the  behavior  of  such  diodes,  without  collisions,  and  with 
few  collisions  f  1  ]  have  been  reported.  These  calculations  as¬ 
sume  that  the  carrier  effective  mass  is  independent  of  energy. 

In  a  device  that  operates  in  a  collision-free,  or  nearly  collision- 
free  regime,  the  carriers  that  contribute  to  the  device  current 
reach  energies  corresponding  to  the  voltages  applied  to  the  de¬ 
vice.  In  small  and  fast  devices,  these  are  of  the  order  of  1  V. 
In  most  semiconductors,  especially  those  with  small  effective 
mass  carriers,  the  bands  depart  from  parabolicity  at  energies 
well  below  1  eV.  The  effective  mass  generally  increases  with 
energy.  Thus  the  current  is  smaller  than  one  would  predict  for 
a  parabolic  band  structure. 

To  assess  the  magnitude  of  this  effect,  we  use  the  conduction 
band  structure  of  a  model  narrow-gap  semiconductor  as  com¬ 
puted  by  two-band  k  ■  p  perturbation  theory.  Such  a  model 
can  provide  a  fair  approximation  to  the  real  band  structure  in 
the  vicinity  of  the  conduction  band  extremum. 

Assuming  a  direct-gap  semiconductor,  and  placing  the  zero 
of  energy  at  the  conduction  band  minimum,  the  theory  gives 
the  hyperbolic  band  structure 


1 

/  2  h2k2 

e  =  —  e. 

1/1  + - 1 

2  * 

The  theory  also  provides  a  relation  between  the  effective 
mass  m *  and  the  bandgap  eg,  but  for  the  present  purpose  these 
two  parameters  may  be  chosen  independently,  for  instance  to 
fit  a  measured  band  structure. 

To  obtain  the  carrier  velocity  as  a  function  of  energy,  we  use 


Collisionless  Space-Charge  Limited  Currents  in 
Semiconductors  with  Nonparabolic  Bands 


3e  _f£k_  I  +  2 
3 k  m*  \ 

and  then  solve  ( 1 )  for  k2 
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Fig.  1.  GaAs  band  structure  with  parabolic  and  hyperbolic  conduction 
band  approximation. 


where  we  have  defined  a  band-structure  limited  velocity 

v,=(etl2m*)lli.  (5) 

Fig.  1  shows  the  current  best  plot  of  the  GaAs  band  structure 
(2|  with  both  (1)  using  tne  experimental  e.  and  m*,  and  a 
/2m*  parabola  superposed  on  the  plot.  The  hyperbola  of 
( 1 )  is  evidently  a  much  better  fit. 

The  Poisson  equation  for  a  semiconductor  uniformly  doped 
with  n0  donors/cm3  carrying  an  electron  current  density  J  is 

ll£.£..22£  +  £  (6) 

rfx  e  e  ev 


where  e  is  the  dielectric  constant  and  u  is  the  electron  velocity. 
Electrons  injected  at  jt  =  0  acquire  a  velocity 


J.  1*1 

*  3k 


(7) 


when  their  energy  €  =  q  V  =  0. 

Using  the  band  structure  gradient  from  (4)  brings  the  Poisson 
equation  into  the  form 


„  _  Q1n0  Jg  (m*vj  -*-0) 

e  eu,  (2 m*u/0  +  02),/J 


(8) 


where  a  prime  denotes  djdx.  Equation  (8)  can  be  integrated 
once  to  give 


(0’)1  -  (0’(O))2  «—  (0*  *  2m*uJJ0),/J  -  ^-^0.  <9) 

ev,  e 

Assuming  space  charge  limitation  0’(O)  =  0,  and  changing  to 
reduced  variables 


s  0  a  J  _ 

2m*  vj  ’  '  qn0v,’  (em*ln0q1)>,i  v. 


(10) 


we  can  write  the  diode  equation,  the  integral  of  (9),  in  the 
form 


dw|>(ir2  +  u)l,i  -  u )  _1/a 


*1 


where  u  ■  v(V  *  F0,  the  applied  voltage)  at  X  *  /(jt  *  L,  the 


Fig.  2.  Space-charge  limited  current-voltage  relation  of  a  collisionless 
semiconductor. 


diode  length)  The  normalized  variables  are:  diode  voltage  in 
units  of  2m*vj  /q -four  times  the  kinetic  energy  the  carriers 
would  have  if  their  effective  mass  were  constant  at  the  band 
edge  value^and  they  were  moving  with  the  limiting  velocity  u,\ 
if  the  k  ■  p  perturbation  theory  is  taken  literally,  this  unit  is 
equal  to  e./q;  current  density  in  units  of  qn0u,-the  current 
that  would  flow  if  the  mobile  carriers  from  the  doping  were 
moving  with  velocity  v,  \  diode  length  in  units  of  i/,rp,  with  Tp 
the  plasma  oscillation  period  of  the  semiconductor. 

Curves  of  current  density  versus  voltage  are  shown  in  Fig.  2 
as  full  curves.  For  comparison  the  current  density-voltage 
curves  one  would  compute  assuming  constant  effective  mass 
are  shown  dashed. 

The  vacuum  diode  analog  3/2  power  law  holds  only  for 
V0  <  eg,  and  then  only  for  1$  1.  That  means  very  short 
diodes— for  GaAs  with  r»o  =  1016,  shorter  than  0.25  pm.  The 
limiting  higher  voltage  behavior  is  J  a  V  rather  than  J  a  V  3/J  ; 
the  J-V  curves  bend  over. 

In  the  typical  size  and  operating  voltage  range  of  GaAs  de¬ 
vices  the  corrections  indicated  here  appear  significant  but  not 
dramatic.  They  are  likely  to  be  rather  more  important  in  de¬ 
vices  made  from  smaller  bandgap,  lower  m*  material. 

We  have  not  extended  the  diode  characteristics  into  the  low- 
current  region  at  the  lower  ends  of  the  curves.  The  problems 
encountered  by  modeling  this  operating  regime  are  rather  com¬ 
plicated,  involving  multiple-valued  functions  and  possibly  spa¬ 
tial  and  temporal  instabilities  [3].  Moreover,  the  neglect  of 
collisions  and  thermal  energies  is  less  justifiable  in  this  low- 
voltage  region. 
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